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Abstract 

We investigate isotropic-isotropic, isotropic-nematic and nematic-nematic phase coexistence in 
binary mixtures of circular platelets with vanishing thickness, continuous rotational degrees of 
freedom and radial size ratios A up to 5. A fundamental measure density functional theory, pre- 
viously used for the one-component model, is proposed and results are compared against those 
from Onsager theory as a benchmark. For A < 1.7 the system displays isotropic-nematic phase 
coexistence with a widening of the biphasic region for increasing values of A. For size ratios A > 2, 
we find demixing into two nematic states becomes stable and an isotropic-nematic-nematic triple 
point can occur. Fundamental measure theory gives a smaller isotropic-nematic biphasic region 
than Onsager theory and locates the transition at lower densities. Furthermore, nematic-nematic 
demixing occurs over a larger range of compositions at a given value of A than found in Onsager 
theory. Both theories predict the same topologies of the phase diagrams. The partial nematic 
order parameters vary strongly with composition and indicate that the larger particles are more 
strongly ordered than the smaller particles. 
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I. INTRODUCTION 



There is a wide range of colloidal particles with platelet-like shape, including materials 
such as gibbsite 1] and certain clays including montmorillonite, laponite and hydrotalcite 
2j-[7|. Clays are some of the most abundant minerals on the Earth's surface and are used 
as pharmaceuticals, cosmetics and catalysts. There is much current interest in the use of 
platelets in nanocomposite materials, e.g. the nematic phase of sterically stabilised gibbsite 
platelets may be used as a template for gibbsite-polymer nanocomposites with nematic 
order Interest in platelet dispersions is also present in geophysics jjj, bio medicine [uj] 
and liquid crystal display (LCD) technology [11 ]. 

Understanding the liquid crystalline phase behaviour of systems of non-spherical particles 



12 



- |l4| is an important topic in modern condensed matter physics. One of the most cele- 
brated cases of a phase transition in such systems is the isotropic-nematic {I-N) transition. 
For athermal model systems, where the particle interactions are hard core, phase transi- 
tions arise purely from entropic contributions to the free energy and the phase behaviour is 
governed solely by density and is independent of temperature. Such models can be used to 
describe lyotropic liquid crystals and phase transitions such as the I-N transition. 

Onsager showed how the formation of liquid crystalline phases can be understood on the 



basis of pair interactions between the constituents of the material [15Hl7| . He considered the 
hard platelet fluid but we know that unlike the case of rod-like particles, his second-virial 
theory does not produce quantitatively correct results for the equation of state and the I-N 
coexistence densities. Onsager himself noted that higher virial contributions are important 
for obtaining reliable results, estimating the ratio B^jB\ at 0(1), with B 2 and B 3 being the 
second and third virial coefficients, respectively. 

Nevertheless, the second-virial theory has been employed to investigate the monodisperse 



platelet system [18|, [19[ . In Ref . 



18] a numerical approach was used to calculate the phase 



diagram of platelets for varying thickness, including the case of zero thickness. This is 
complemented by a calculation for the equation of state in both the / and N phases for van- 
ishingly thin platelets in Ref. 19|. The first off-lattice simulation study of the I-N transition 
for infinitely thin platelets was carried out in Ref. j^. This showed that the phase diagram 
differs substantially from the Onsager prediction and that the I-N transition is actually 
much more weakly first order and occurs at lower densities than predicted theoretically. 
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The authors also carried out a fifth order virial calculation for the equation of state. More 



accurate predictions of the higher virial contributions for disks were presented in Ref. 21], 



where simulation results are reported for hard cut spheres and more rece ntly in Ref. 
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Later simulation work was carried out on polydisperse platelet systems [23] and systems 
of platelets with different polygonal shapes (e.g. hexagons, triangles) 24j. Further simula- 
tion results of model circular platelets were reported in Ref. [25] and simulations alongside 
an integral equation approach for mixtures of rods and disks were carried out in Ref. j^ . 
Simulations of binary platelet systems have not yet been carried out. 

Binary mixtures of particles of different shape and/or size are interesting due to the 
richness of the phase diagrams they may exhibit. Binary rod mixtures form a prominent 
example. Studies include mixtures of thick and thin rods [27] and long and short rods 
using Onsager theory as well as using Parsons-Lee scaling |29|. The phase behaviour 



in binary mixtures can include: the fractionation effect, whereby the larger particles go 
preferentially into the nematic phase; widening of the biphasic region; a re-entrant / — > N — » 
/ phenomenon on increasing density; the possibility of demixing into two different isotropic 
states and/or two different nematic states and triphasic equilibria (see e.g. Ref. H for 
examples of these phenomena). Nematic- nematic (N-N) demixing, at high enough pressures, 
can be viewed as a result of competition between orientational entropy of the smaller platelets 
favouring mixing, and the entropy of mixing j^. The N-N phase separation for binary 
mixtures of rods, including the high density regime, is studied in detail in Ref. 30]. 



Binary mixtures of thin and thick platelets have been investigated 3lJ, |65] with the 
Parsons-Lee scaling of the Onsager functional 32j-|34|. Studies based on the Zwanzig model 
for binary hard platelets, where the particles are restricted to occupy only three mutually 
perpendicular directions, have been carried out for the bulk and interfacial properties of the 
demixed phases. Rich phase diagrams, involving isotropic and nematic phases have been 
reported in Refs. 351-137]. A recent review 38] of platelet fluids contains a summary of these 
results. Recently Verhoeff et al. 39] have investigated experimentally and theoretically the 
phase behaviour of colloidal platelets with bimodal shape distribution. Their theory is based 
on the Onsager- Parsons free energy and a cell approach for the columnar {Col) state 40]. 
The authors find agreement between their experimental findings and theoretical predictions 
for sufficiently large thickness ratios. The phase diagram features an I-N density inversion 
and triphasic I-N -Col equilibrium. 



3 



Fundamental measure theory (FMT) is an approximate non-perturbative density func- 
tional theory (DFT) 41], originally proposed by Rosenfeld for additive har d sp here mixtures 
The approach was later generalised to other convex shapes [44|, |45| , which led to 
subsequent work Q, 4?|. The bulk I-N coexistence densities (scaled by the cube of the 
platelet radius) and nematic order parameter at the transition (cj, cjv and Sn, respectively) 
were calculated by Frenkel and Eppenga in Ref. [J by simulation; for more recent sim- 
ulation results see Ref. 481. The values previously obtained from FMT are cj = 0.418, 
Cn _ = 0.46 and Sn = 0.53 [48] . which are in agreement with the present study. Recently 



49] these values were improved using the same method but with increased resolution to 
cj = 0.419, cn = 0.469 and Sn = 0.533 [50j. The FMT functional for pure platelets was 



later utilisec 
a hard wall 



to study inhomogeneous situations including the I-N interface and wetting at 
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51] and capillary nematisation of platelets between two parallel walls 48]. 



Generalising the theory for the corresponding one-component system 52] , we here propose 
a functional to describe binary mixtures of vanishingly thin circular platelets. Our theory 
features the exact virial second order term in density and an approximate term of third order 
in density. We investigate three types of demixing phase behaviour in the case of binary 
platelets with varying size ratio, finding I-N and N-N phase coexistence. We do not find 
stable I- 1 demixing (as could be driven by the depletion effect [53] ) for the regimes considered 
in the present work. We restrict our attention to uniaxial arrangements of the (uniaxial) 
platelets, as we do not expect biaxial arrangement of the particles to occur. We study a 
range of size ratios in this investigation, ranging between A = 1.1 and 5. We present the 
phase diagrams in different representations to facilitate comparison with simulations which 
may be performed in different ensembles or experiments. We expect the phase diagrams 
from FMT to be quantitatively more accurate than those from Onsager theory, which we 
calculate as a reference. The topologies of the phase diagrams are the same in both theories 
for our chosen values of the size ratio between two species. Although the integral kernel, 
which represents the pair excluded volume term, is the same for long thin rods as it is 
for platelets of vanishing thickness, the results from Onsager theory cannot be obtained by 
simple scaling of literature results for binary mixtures of rods. 

This paper is organised as follows. In Sec. |TT] we define the model, outline the density 
functional theory and the conditions for thermodynamic stability and phase coexistence. In 
Sec. II III we present results for the phase behaviour of binary platelet mixtures. In Sec. IIVI 
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we provide conclusions and an outlook on possible future work. 



II. DENSITY FUNCTIONAL THEORY FOR BINARY HARD PLATELET MIX- 
TURES 

A. Pair Interactions and Model Parameters 

We consider a binary mixture of hard circular platelets with vanishing thickness and 
continuous positional and orientational degrees of freedom. Species 1 and species 2 have 
radii Ri and i? 2 respectively, with R2 > R\. The pair potential Uy between two particles i 
and j, where i, j = 1, 2, is infinite if the geometrical shapes of the two platelets overlap and 
zero otherwise and is hence given by 



(00 if particles overlap 
(!) 
otherwise, 

where r and r' are the positions of the particle centres and u) and a/ are unit vectors 
indicating the particle orientations (normal to the particle surface). The size ratio 



characterises the radial bidispersity and is the only control parameter in the model. We 
characterise the thermodynamic state by two dimensionless densities C\ = p\R\ and c 2 = 
P2-R1, where pi and p 2 are the number densities of the two species, p« = Ni/V, where iVj is 
the number of particles of species i — 1,2 and V is the system volume. The composition 
(mole fraction) of the (larger) species 2 is x = p2li.P1 + P2) and the total dimensionless 
concentration is c = R\ (pi + p 2 ) = C\ + c 2 . 

B. Grand Potential Functional and Minimisation Principle 

Density functional theory (DFT) is formulated on the one-body level of the density dis- 
tributions pi(r,u:). The variational principle |4l| states that the true equilibrium density 
profile is the one which minimises the grand potential functional Q and so obeys 




(2) 



(5^([pi,p 2 ],pi,p 2 , V,T) 
Spi(r,u) 







(3) 
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where pi is the chemical potential of species i = 1,2 and T is absolute temperature. The 
grand potential functional can be decomposed as 



fi([pi,p 2 ],/ii,/i2, V,T) = F id ([p h p 2 ],V,T) 

2 

+ F exc (iPup2],V,T) + J2 JdrJ dw/j^Cr.w)-^), (4) 



8=1 



where the spatial integral (over r) is over the system volume V and the angular integral 
(over uj) is over the unit sphere; V^ xt (r, oj) is an external potential acting on species i\ 
F e xc([pi, P2], V, T) is the excess (over ideal gas) contribution to the Helmholtz free energy 
and describes the inter-particle interactions. The free energy functional for a binary ideal 
gas of uniaxial rotators is given by 

2 

0F id {[ Pl ,p 2 ],V,T) =J2 [ dr [ dwA(r,«) 

i=i •* •* 

x [ln( A (r,a;)Af)-l], (5) 

where Aj is the (irrelevant) thermal wavelength of species % and /3 = l/(fcsT), where ks 
is the Boltzmann constant. We let Aj = Ri, which is equivalent to fixing an additive 
constant to the chemical potential and hence does not affect any observable properties of 
the system. One systematic way to express the excess free energy functional is to expand it 
in a virial series in density Onsager theory is based entirely on the second- virial level. 
FMT (as described in Sec. Ill Dl below) approximates higher order terms using single particle 
geometries. Nevertheless we find it useful to give the terms in the virial expansion up to 
third order in density explicitly: the second and third order contributions to the (exact) 
virial series for the excess free energy (3F exc ([pi, p 2 ], V,T) are given respectively by 



| [I + I +=*f]. (6) 



1 2 J 
1 2 2 



4[L + L+ 3 I\+ 3 I\]- < 7 > 



'"11 22 11 12" 



where each line in the diagrams represents a Mayer function /y(r — r', u>, a/) = 
exp (— f3iiij(r — r', a?, uj')) — 1, which equals —1 if the two particles overlap and zero oth- 
erwise. The shaded circles, field points, indicate multiplication by the one-body density 
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Pi(r,Lj) and integration over the coordinates r and uj [16]. The number alongside each field 
point represents the species i. Here we consider only spatially homogeneous states, such that 
the one-body density does not depend on position r and, for the case of uniaxial nematic 
states considered in this paper, depends only on the polar angle 9 of u> with respect to the 
nematic director. Hence the one-body densities factorise as Pi(r, uj) = p^i{9) where ^i(9) is 
the orientational distribution function (ODF) and pi the number density of species i. Before 
laying out the FMT we focus on the second-virial level. 



C. Onsager Second- Virial Theory for Binary Platelets 

The diagrams in Eq. (Q, using the fact that pi(r, uj) = pi^i{9) for spatially homogeneous 
states, become 



j 

- II = -\piPij du J du>'j drMwu^VfflVjffl, (8) 

i 

where we have renamed r — r' — > r. The spatial integral over the Mayer bond yields (minus) 
the excluded volume —£ij(uj,uj') between two particles of species i and j, as a function of 
the angle 7 between u) and uj'. Hence 

Eijiyj, uj') = - J dr/u(r, w, uj') = 2n(R 2 i R J + R]^) sin 7 , (9) 



Therefore 



3 r% 

- * = 167r 2 (i2?iL- + RlRjpiPj / 2 d9 sin 9 (10) 

V* ' Jo 

1 



x / d9' sm9'K(9,9 , )^ i (9)^ j (9'). 
Jo 

The integrals in Eq. fflQ|) (omitting the prefactor) can be written as 



2-7T r~K p2n 

d9 sin 9 / dd>' I d9 sin 9' 



JO Jo Jo 

x * 1 (0)* 1 (0 / )sin7, (11) 

where 9 and 9' are the polar angles of two platelets with respect to the nematic director 
and and <$' are the azimuthal angles. Due to the inversion symmetry of the nematic 
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state Jq d9 — 2 j Q 2 dO. In order to deal with the azimuthal integral we introduce the kernel 
K(9, 9') via 



K{9,&) 



2- 



>sm7 



2tt 



i\2 



i-2n 

/ d(fiyl — (cos 9 cos 9' + sin 9 sin 9' cos (f)) 2 , 
Jo 



(12) 



where we have renamed <fi — <fi' — > <p as the difference between the azimuthal angles of the 
two platelets. Adding all three terms and multiplying by Rf yields the Onsager contribution 
to the excess free energy in the dimensionless form 



(2) 



V 



—Rl = 167r 2 c 2 d9sin9 d9' sin 9' K{9, 9' 



(1 - xf^^^ff) + x 2 \ s y 2 {8)y 2 {6') 
+ x(l - x)(A 2 + A)*i(0)tf 2 (0') 



(13) 



where the superscript of Fixe represents the order in density of the excess free energy. 



D. Fundamental Measure Theory for Binary Platelet Mixtures 

We generalise the monodisperse functional to the case of binary mixtures using an 



ap- 



proximate term at third order in density which is based on the FMT developed in Ref. [52 ]. 
Our theory contains the exact second order Onsager term and instead of using any higher 
order terms from the series, such as the exact third virial level (jTj), an approximate term 



which is of third order in density is used 54j. This term is nonvanishing (and constant) for 



cases with common triple intersection of the three platelets involved. There are no higher 



order terms due to the scaled-particle roots 55J of the approach; the vanishing volume of 



the platelets truncates the series. Global prefactors are used to compensate for lost cases 



56l . 1571 ] . We postulate the excess free energy 

PF exc ([ Pl ,p 2 ],V) = J dr J du J dw' 



< D (r,u;)n°(r,a;') 



+ ^n?(r, u,K DD (r, w, u;'K(r, , (14) 

(where the right hand side is independent of T). The first term of the sum in Eq. (TBI) is 
equivalent to the Onsager contribution to the excess free energy and the second is the FMT 
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contribution. The weighted densities are related to the bare one-body densities, pi(r, u>), via 



2 r 

? D (r,o;) = ^ / ^< a (r, W »* ft (r )W '), 
i=i J 

2 
i=l 

= E/ ^V Da (r,"VVO*Pi(rV'), 



DDD / 

n 2 (r, u>,u> 



(15) 
(16) 
(17) 



where * represents the three-dimensional convolution h(r) * g(r) = J d 3 xh(x)g(jx. — r). We 



have kept the notation of Ref. 



52] where the upper index D (disk) is indicative of the number 



of particle orientations that appear in the weight function (below) or weighted density. The 
weight functions for species i are given by 



w™(r,u>) = 6{R i -\r\)6(r-u>)/8, 
w$ i (r,u) = 2e(R i -\T\)5(T-e>), 



,DDi , 



R, 



,DDDi, 



|u> • (u> x r)\w l [r,ut), 



wz~'(t,w,w',u/") = -\uj ■ (w' x u")\w^{t,u), 

71 



(18) 
(19) 
(20) 

(21) 



where S(-) is the Dirac delta distribution, G(-) is the Heaviside step function and x denotes 
the vector product, such that u) ■ (u/ x u>") is the triple scalar product. Note the modulus 
in Eqs. f l20|) and (j2ip . For spatially homogeneous states, Eq. ( I14p becomes 

f3F exc ([pi,p 2 ]) 



V 



duj I du)' 
1 



,DD / 



n 1 \U))n 2 [OJ , 



(22) 



Inserting the definitions of the weighted densities (IT5"l) - (ITTj) into the excess free energy (Tfj 
we obtain 

12 2 

• r r r r \u-(u r xu")\ 

(23) 



/3F exc ([ Pl ,p 2 ]) = -]- [T + T +2T] + fdx f du j du>' f du; 

12 1 J J J J 



3tt 2 



xn°(x, w)n 2 D (x, w')n£(x, w") 
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The fundamental measures of a platelet of species i are the integral mean curvature £ 4 - MC 
TrRi/4 and the surface £f = 2wRl The first term of Eq. ([22} may be expressed as 



A\u),u}') = - 



+ (^ M % s + e% S )^iW^(^)/2 



sin 7, 



(24) 



remembering that ^(u;) = ^i{9) for uniaxial nematics. Eq. ff24l) is simply the fundamental 
measures interpretation 52J of the Onsager contribution to the free energy. Note that the 
term which represents the excluded volume between two particles of different species is given 
by £2 MC £i + ^} MC ^2 • This leads to a scaling of the excluded volume by A 2 + A, which is quite 
different from the scaling that occurs for binary mixtures of rods. The second term in 
Eq. (I22p is the FMT contribution to the excess free energy. This is given by 



(3) 
exc 



V 



du / dJ \ dJ'B'iu^u'J 1 



(25) 



where 



B\u>,u>',u)") 



|u? • (u)' x uj") 

3^ 



(26) 



We choose coordinates such that uj = (sin 9, 0, cos 9), uj' = (cos 9' sin 9', sin <ft' sin 9', cos 9') 
and uj" = (cos 9" sin 9", sin 0" sin 9", cos 9") where 9,9' and 9" are the polar angles of three 
platelets. The third order contribution in density to the FMT excess free energy in these 
coordinates is given by 



p(3) 108-2 /-f /•§ /•§ 

rexc ^3 = i^ c 3 / d9sinQ de > sin6 > d9" sin 9" L(9, 9', 9") 
V 3 y J J 

[(1 - x) 3 ^)^')^!^") + X 3 \^ 2 (9)^ 2 (9')^ 2 (9") 

+ 3x(l - x) 2 \ 2 ^ 1 (9)^ 1 (9')^ 2 (9") + 3x 2 (l - x)X 4 ^ 1 (9)^ 2 (9')^ 2 (9")], 

(27) 
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where the kernel L(9, 9', 6") is 

r-2n p2n 



L(6, 9', 0") = / / d(//d<f>"\u ■ (u' x u 
Jo Jo 

"2-k p2n 

j it i ill 



sin #(sin 0' sin 9' cos 9" 

o Jo 

+ cos 9' sin 0" sin 9") + cos #(cos 0' sin 9' sin 0" sin 9" 
-sin0'sinfl'cos0"sinfl")|. (28) 

The full form of the excess free energy, as used in the calculations described below, is given 
by the sum of Eqs. f|T3|) and (J2"T|) . In practice, along with Eqs. (}T2"]) and (1251) . these require 
numerical computation on a grid as described in the following. 

E. Self-consistency equations for the orientational distribution functions 

The minimisation principle ([3]) together with the FMT approximation (122]) for 
F exc ([px, p 2 ]) leads to two coupled Euler-Lagrange equations for the ODFs: 



# 2 (0) = ^exp 

^2 



7T 

- 8ttc / 2 d0' sm9'K(9, 9')[(1 - + \x{\ 2 + A)* 2 (0')] 

7T 7T 

-32ttc 2 [ 2 d9'sm9' d9" sin 9" L(9, 9' ,9") 
Jo Jo 

x [(1 - xf^> 1 {9')^ 1 {9") + 2x(l - x)A 2 ^ 1 (0 / )^ 2 (O + x 2 \^ 2 (9')^ 2 (9")} 

(29) 

rf 1 
-8ttc / ^ , sin^(^,^)[^A 3 ^ 2 (^) + 7j(l-a;)(A 2 + A)^i(^)] 

-32ttc 2 J* d9'sm9' j d9"sin9"L(9, 9', 9") 
Jo Jo 

x [x 2 A 6 ^ 2 (^)^ 2 (^) + 2x(l - x)\^ 1 (9')V 2 (9") + (1 - s) 2 A 2 *i(0')*i(0")] 

(30) 

where the constants Z\ and Z 2 are such that the normalisation j doj^^uS) = 1, for i — 1, 2. 
Neglecting terms of order c 2 in the exponentials, the equations for FMT reduce to the 
Euler-Lagrange equations of Onsager theory. We solve Eqs. ( 129|) and (1301) numerically with 
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a straightforward extension to the iterati ve p rocedure given in Ref. 58] and numerical tech 



niques similar to those described in Ref. 59]. The 9- and 0-grids are defined on [0, 7r/2] and 
[0, 27r] respectively. The #-grid is divided into 200 equal steps on [0,7r/8], where the ODF 
changes most rapidly in the nematic phase and into 50 equal steps on [7r/8,7r/2] where the 
ODF is almost zero. The 0-grid is divided into 200 equally spaced intervals on [0,27r]. We 
start the iteration with two initial trial distributions, for example a normalized Gaussian 
(c/7r 2 ) exp [— 2c 2 8 2 /ii] or a constant distribution 1/4-7T. The choice of the constant distribu- 
tion is more efficient at low densities where the system is expected to be isotropic. These 
guesses are substituted into the right hand sides of Eqs. ( 1291) and ( |30l) to obtain a new pair 
of ODFs, \I/j inew (6!). This procedure is repeated until max \^i tnew (0) — ^i )0 id(^)| < t, i = 1, 2, 
where t is the tolerance given by the magnitude of the largest acceptable deviation to the 
ODF of the previous step: \E f i, id(^)- We set t = 10~ 9 . For FMT the solutions take a longer 
time to converge than for Onsager theory, as each step has a slower execution time due to 
the increased complexity of the coupled equations (l2~9l) and (l30~j) . 

F. Conditions for Phase Coexistence 

Once we have found the ODFs we solve the phase coexistence equations. The require- 
ments for phase coexistence between two phases A and B are the mechanical and chemical 
equilibria between the phases as well as the equality of temperature in the two coexisting 
phases (which is trivial in hard-body systems). Hence we have the non-trivial conditions 

P A = P B (31) 

and 

nt = if, (32) 

where i — 1, 2 again labels the species. We calculate the the total Helmholtz free energy F 
numerically by inserting ^i{9) into Eqs. (jHJ), (fl3l) and (|27j) . Likewise, the pressure can be 
obtained numerically as 

F , * _ d(F/V) 
i=i 



P=-7+E^ (33) 



and the chemical potentials as 



d{F/V) 

W = -fa-' (34) 
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We define a reduced pressure p* = f3pR\ and reduced chemical potentials //* = Eqs. (I3T]) 
and (1321) are three equations for four unknowns (two state points each characterised by two 
densities). Therefore, regions of two-phase coexistence depend parametrically on one free 



parameter. Eqs. (13"3"1) and (13~4|) are solved numerically with a Newton- Raphson procedure 60]. 
The resulting set of solutions yields the binodal. I-N-N triple points are located where the 
I-N and N-N coexistence curves cross. Therefore, regions of two-phase coexistence depend 
parametrically on one free parameter (which can be chosen arbitrarily, e.g. as the value of 
concentration x in one of the phases). 



G. Equation of State for the Isotropic Phase 



Analytic expressions for the free energy, pressure and chemical potentials for the isotropic 
phase may be found on insertion of ^(0) = 1/(4-71") into the ideal ([5]) and excess fl22|) parts 
of free energy functional. For this purpose, we use 

-7-2 
Xl 7T 



n Jo 



57T 



and 



d6sm6 / dO'sinO' / d6" sin 9"L(6, 6', 9" 
1 1 J o Jo 



Z2_^ 
167T 2 



(35) 



(36) 



where the integrals X\ and Z2 are calculated in Appendix A. The expressions for the FMT 
isotropic free energy, (3F iso /V ', pressure, p* so and chemical potentials /x* iso are 



V 



c(l — x) 



mo Rl = c(l - x) In ( ~' \ + cxln ( ^- ) - c 



47T 



CX \ 

4n~) 



+ yc 2 [(1 - xf + x 2 A 3 + x{\ - x){\ 2 + A)] 



3 



(37) 



C 2 7T 2 



c + [(1 - s )2 + X 2 A 3 + (A 2 + A)(a , _ ^ 



+ 



2vr 2 c 3 



z 3 A b + 3(x 2 - x 3 )A 4 



+ 3(x - 2x 2 + x 3 )X 2 + (1 - x 



(38) 



Ki so = In (^^) + IT t 2(1 ~ X) + X(A ' + A) ] 
+ cV[(l - x) 2 + 2z(l - a;)A 2 + x 2 A 4 ], 



(39) 
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^iso = In (^) + ^ [2xA 3 + (1 - x)(A 2 + A)] 

+ c 2 tt 2 [x 2 X 6 + 2x(l - x)A 4 + (1 - xf\ 2 }. (40) 

The Onsager versions of these equations are given by the same expressions but without the 
final bracketed term in c 3 for f3F- lso /V and p* so and c 2 for //* iso . The authors of Ref. [61 ] 



claim that inclusion of the exact third virial term along with the Onsager term, at least for 
monodisperse platelets, would give a worse equation of state in the isotropic phase than the 
Onsager term alone. 



H. Isotropic-Nematic Bifurcation Analysis 

On increasing the density in the isotropic state, a point is reached known as the bifurcation 
density, where an infinitesimal nematic perturbation destabilises the system. The first I-N 



bifurcation analysis for a liquid crystalline system was performed in Ref. 



. 3. 



This was 



extended to a class of liquid crystal models in Ref. [63]. The bifurcation concentration lies 
inside the coexistence region for the monodisperse case of platelets [26] but as we will see, 
this is not always true in the binary case. We insert *&i(9) = [l + eiP2(cos#)]/47r and ^2(6) — 
[1 + e2-P 2 (cos #)]/47r into the free energy ([5]) and (122]) . where €\ and e 2 are small parameters 
measuring the strengths of the nematic perturbation and P2(cos8) = (3 cos 2 (#) — l)/2 is the 
second degree Legendre polynomial in cos#. We then extract all second order terms, i.e. 
those proportional to e\ , el and €±€2, respectively. The coefficients of these terms are denoted 
by cii(c), a 2 (c) and ai 2 (c), respectively. These involve integrals of Legendre polynomials and 
are obtained using the integrals X 3 ,X 4 ,X 5 defined in Appendix A. One then solves det M = 



53j where 



Here 



M=| ai(c) " 12(C)/2 ]. (41) 
a i2 (c)/2 a 2 (c) 



a 1 (c)=^(l-x)-^c 2 (l-x) 2 

-^c\l-xf-^x(l-x) 2 \ 2 , (42) 
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7f2 A 6 c 3 a; 3_^! c 3 ;r 2 (1 _ a;)A 4 ) (43) 



40 40 

ai2(c) 7T 2 2 2 

-c x(l - x)(X + X) 
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—c 

Id' ' 10 



_ ^! c 3 x(l _ X )2 A 2 _ ^! C 3 X 2 (1 _ ^4 (44) 



The results for the spinodals were checked by running the self-consistency program for the 
solutions of the coupled ODFs (129]) and (130]) with the trial functions l/47r; the locus of c 
values (for a given value of A) for which the maximum number of iterations occured was 
found to agree numerically very well to the I-N spinodal. To calculate the spinodals in the 
(x,p*) representation we insert the (x, c) values which form the spinodal into the expression 
for the isotropic pressure (|38p . In the monodisperse limit, one needs to solve the following 
cubic polynomial to calculate the bifurcation point c* 

TT 2 9 7T 2 , , , 

c * - Y c * ~ T c * = °' (45) 

where the bifurcation concentration is c* = 8/tt 2 = 0.811 for Onsager theory and c* = 0.434 
for FMT, the latter much closer to the limit of stability of the isotropic phase as observed 



in simulations 



26) 



I. Symmetry-Conserved Demixing Spinodals 

A thermodynamic phase is locally stable if the determinant of the Hessian matrix of the 
Helmholtz free energy density with respect to the species densities is positive. The spinodal 
is the limit of stability, defined by det N = 0, where 

( d 2 (F/V) d 2 (F/V) \ 
d 2 (F/V) d 2 (F/V) J ' 1 ; 

dpidp-2 dp2 / 

For the 1-1 spinodals, we insert ^i^) = ^2^&) = 1/47T into the free energy ( 1371) . This yields 
an analytic solution of Eq. (146]) for Onsager theory, given by 



1\ / 1 \ vr 4 
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and for FMT by 

(V 2 + 27r 2 pi + 2p 2 n 2 X 2 + — 
V " pi 

x (Va 3 + 2p 2 7r 2 A 6 + 2 Pl vr 2 A 4 + - 



P2- 

(2 Pl 7r 2 A 2 + 2p 2 7r 2 A 4 ) 2 = 0. (48) 



I-I demixing never occurs for the values of A considered here; indeed solutions of Eq, 
only begin to exist at about A = 20. For the N-N spinodals, there are no such analytic 
equations (for the Zwanzig model, see Ref. 3^]). For practical reasons we rather solve 

dp* 



dp2 



0, (49) 



/'i 



which is equivalent [53| to solving det N = 0. In order to calculate the N-N spinodals we 
have to numerically evaluate the left hand side of Eq. ( I4"9l) . Exchanging the species labels in 
Eq. (14"9"|) one obtains the same results. The N-N spinodals are calculated to give an idea of 
the location of the N-N phase boundaries. For cases where there is N-N coexistence closed 
by critical point, the spinodal and the binodal coincide at the critical point. 

III. RESULTS 

We have calculated the phase diagrams of binary platelet mixtures for seven different 
size ratios, A = 1.1, 1.4, 1.7,2,2.5,4 and 5. Fig. la shows the results for the slightly asym- 
metric case A = 1.1 in the (x, c) representation. The I-N transition is the only type of 
transition which we find for this size ratio. The FMT results show that the I-N transition 
concentrations at x = are cj = 0.418 and cn = 0.46 respectively, in agreement with the 



monodisperse results [25]. At x = 1, the coexistence values are cj = 0.418/A 3 = 0.314 
and cjy = 0.46/A 3 = 0.346. The coexistence curves interpolate smoothly from x = to 
x = 1. The tie lines joining coexisting isotropic and nematic phases are naturally verti- 
cal at x = 0, 1 (corresponding to the cases of the pure systems of small and big platelets, 
respectively) whereas between x = and x = 1 they vary in gradient, leaning with large 
positive gradient as x increases from zero composition, less so at about 50% composition 
and then more so again on approaching x — 1. Therefore, there is stronger fractionation at 
50% composition than towards x = and x = 1, such that the isotropic phase is dominated 
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by particles of species 1 (smaller species) and the nematic phase by particles of species 2 
(larger species). The biphasic region, where there is coexistence between the isotropic and 
nematic phase, is very narrow in FMT. At zero composition, Onsager theory predicts a 
density jump of ~ 22% whereas this is only ~ 9% for FMT, which agrees more closely to 
simulation results of 8% 61J. The bifurcation concentration for FMT is c* = 0.434 at x = 
and c* = 0.434/A 3 = 0.326 at x = 1. The spinodal for FMT lies closer to the isotropic 
phase boundary than the nematic phase boundary (whereas the converse is true for Onsager 
theory). Fig. lb shows the same results but in the (ci,c 2 ) representation. The results for 
both theories again interpolate smoothly between the two pure limits. The two branches 
of the binodal in this representation move from the Ci-axis to the C2-axis with increasing 
composition. Hence the tie lines move from being horizontal on the Ci-axis to vertical on 
the C2-axis. In the (x,p*) representation (Fig. lc) the tie lines are horizontal due to the 
requirement of equal pressure in the coexisting phases. The spinodal lies above the nematic 
branch of the binodal in this representation because it is obtained by inserting the bifur- 
cation densities into the isotropic equation of state, which yields a higher pressure than 
the coexistence value. The binodal obtained from FMT is located at significantly smaller 
densities as compared to that from Onsager theory. The isotropic end of the tie line is 
at a lower composition than the nematic end of the tie line. For Onsager theory at zero 
composition, the isotropic branch of the binodal intersects the c-axis at Cj = 0.666 and the 
nematic branch intersects at cn = 0.849, in agreement with the monodisperse limit found in 
earlier work [52]. The binodals interpolate smoothly from x = to x = 1 where the values 
of c corresponding to the I-N coexistence concentrations are cj = 0.5 and cn = 0.638. The 
nature of the tie lines is similar to FMT. The I-N spinodal lies between the two branches 
of the binodal and interpolates smoothly from zero composition, where c = 8/n 2 = 0.811 to 
x = 1 where c* = 0.609. 

The nematic phase of a mixture of two components can be characterised by two partial 
nematic order parameters, Si and S2 defined by 

Si = 4tt [ 2 d0sin0^(0)P 2 (cos0). (50) 
Jo 

The total nematic order parameter is the weighted average 

S tot = (l-x)S l + xS 2 . (51) 
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Here we investigate the behaviour of these quantities at I-N and N-N coexistence. We have 
chosen S to t to be a simple weighted average of the Si. As such, each particle contributes 
to St t independently of its size. Of course there are other suitable choices which may be 
appropriate to certain applications; the Sj in the sum for the total order may be weighted 
by the surface area of the platelet, for example. Due to the absence of particles of species 
2 at x = 0, Stot and Si take the same value. Similarly at x = 1, S to t and S2 take the same 
value. As x increases from to 1, S to t changes smoothly. In Fig. Id the FMT values for Si, 
S2 and Stot are smaller than those obtained from Onsager theory, which places the values 
much higher. However FMT predicts that the difference between Si and S2 to be bigger 
than Onsager theory does. S to t at x = takes the same value S to t = 0.531 as previous work 
(and agrees well with simulation 25]) in the monodisperse case and the the order parameters 
vary smoothly in the same manner as for Onsager theory. 

In Fig. 2 we show the results for A = 1.4. There is a widening of the biphasic region 
between the two pure components, clearly seen in Fig. 2a where we show the phase diagram 
in the (x, c) representation. FMT again predicts that the mixture undergoes I-N phase 
separation at lower densities than Onsager theory. The tie lines become less steep with a 
smaller positive gradient than for A = 1.1 for intermediate values of composition indicating 
that there is larger difference in mole fraction between coexisting isotropic and nematic 
states. The widening of the biphasic gap is even more noticeable in the (ci, C2) representation 
(Fig. 2b), especially for the case of Onsager theory. In the (x,p*) representation (Fig. 2c), 
the spinodal again lies above the nematic binodal for Onsager theory. However, for FMT 
we see that while the spinodal is above the nematic branch of the binodal close to x = 
and x = 1 the curve enters the biphasic region in between about x = 0.1 and x = 0.6. 
Si is smaller than S2 and their difference has increased for both theories, with FMT still 
possessing the larger difference suggesting that the particles of species 2 are significantly 
more ordered than species 1 for a given mole fraction at coexistence. 

In Fig. 3 we plot the graphs for A = 1.7. The I-N biphasic gap becomes even more 
pronounced. The results for FMT (Fig. 3a) show that tracing along the nematic branch 
of the binodal as x increases leads to an increase in c. At approximately x = 0.2, c then 
decreases and near x = 0.6 bends back on itself before approaching x = 1. This bending of 
the binodal constitutes a re-entrant phenomenon. There is also a large range of compositions 
between about x = 0.1 and x = 0.6 for which the biphasic I-N phase overlaps for both 
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theories, which is also seen clearly in the (ci,c 2 ) representation (Fig. 3b). In the (x,p*) 
representation (Fig. 3c) both theories predict a region of composition values for which the 
spinodal lies inside the biphasic region. For illustration of the re-entrant part of the phase 
diagram one should keep a constant fluid composition near x = 0.6 but increases the pressure 
from p* = to p* = 1: the state changes from / — > I+N 2 — > N — > I+N 2 — > N where N 2 is 
a nematic phase composed mostly of particles of species 2. The order parameters along the 
nematic branch of the binodal are shown in Fig. 3d. Si remains lower than S 2 for a given 
composition, again suggesting that the particles of species 2 are more ordered than those 
of species 1 at coexistence. When the re-entrant feature occurs at about x = 0.6, Si, S 2 
and S'tot drop sharply. The difference between Si and S 2 along the nematic branch of the 
binodal has become much larger than for A = 1.4. 

In Fig. 4 we present the phase diagrams for A = 2. The most striking feature after 
increasing the size ratio to A = 2 is the stable N-N coexistence region. For FMT the (ci, c 2 ) 
representation (Fig. 4a) reveals that this region is in the form of an almost symmetric hump 
suggesting that the fractionation between the two distinct nematic phases becomes less with 
increasing density. There is also a triple point which is in the form of a triangle connecting 
a low composition isotropic phase, a nematic phase composed mostly of particles of species 
1 (Ni) and a nematic phase composed mostly of particles of species 2 (N 2 ). The symmetric 
N-N coexistence region is more clearly seen in the (x,p*) representation (Fig. 4c) between 
x just greater than and x = 0.45 ending in a critical point at (x,p*) = (0.19,2.72). 
In the (x,p*) representation the triple point collapses onto a line and the pressure at the 
triple point is approximately p* = 1.75. The I-N biphasic region has become increasingly 
pronounced. The topology of the phase diagrams is the same as that obtained in Onsager 
theory results (Fig. 4b and Fig. 4d). The re-entrant feature obtained from Onsager theory 
is less pronounced than in FMT. The N-N phase separation occurs over a larger range of 
composition values (up to near x = 0.5) and the triple point occurs at approximately a 
unit of reduced pressure higher than the triple point predicted by FMT, nevertheless, as we 
emphasise, the topologies predicted by both theories are the same. The I-N spinodal enters 
the biphasic region for both theories. We postpone the discussion of the order parameters 
for A = 2 and higher size ratios until the end of this section. 

In Fig. 5 we show the results for A = 2.5. There is N-N demixing in the system, as is 
observed for A = 2, however there is a big difference in the toplogy of the phase diagrams in 
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that the demixing now does not end in a critical point. For FMT in the (c±, c 2 ) representation 
(Fig. 5a) the N-N region opens up and the two branches of the binodal extend outwards 
suggesting that the demixing extends to arbitrarily high density. No critical point is observed 
up to the densities we examine; we follow the phase boundaries up to c = 1 for FMT and 
c = 1.9 for Onsager theory. A similar splaying of the phase boundaries is observed in Fig. 5b 
in the (c±, c 2 ) for Onsager theory, suggesting that the topology obtained with FMT is correct. 
Also, the I-Ni region has shrunk to a tiny region close to zero composition, making the triple 
region mostly dominated by the nematic state rich in species 2. In the (x,p*) representation 
for FMT (Fig. 5c) the re-entrant feature has become extremely pronounced, extending as 
low as approximately x = 0.5 in the (x,p*) representation and the triple point occurs at 
approximately p* = 1.8. In the (x,p*) representation for Onsager theory (Fig. 5d), the re- 
entrant feature has also become more pronounced but less so than for FMT. The triple line 
extends over approximately the same composition range (x < 0.8) as for FMT but occurs 
at just over a unit of reduced pressure higher than for FMT. 

In Fig. 6 we present the phase diagrams for A = 4. In the (01,02) representation for 
FMT (Fig. 6a) the strength of the I-N2 fractionation effect becomes very large. The N-N 
separation is very wide giving a huge immiscibility gap, hence we do not show it. Fig. 6b also 
shows a large immiscibility gap in the (01,02) representation for Onsager theory. We again 
have confidence that the results from FMT are quantitatively more accurate than those from 
Onsager theory. At concentrations c\, 02 < 0.1 the coexisting compositions already approach 
closely the ci,c 2 axes highlighting that coexisting compositions are close to x = 0, 1 as is 
clearly seen in the (x,p*) representation (Fig. 6c for FMT and Fig. 6d for Onsager theory). 
The general trend of the re-entrant bend moving to lower composition with increasing A has 
continued here, reaching as low as about x = 0.3 for FMT and near x = 0.5 for Onsager 
theory. 

In Fig. 7 we plot the results for A = 5. In the (ci,c 2 ) representation for FMT (Fig. 7a) 
the I-N2 coexistence region is very pronounced. The larger platelets have an area twenty- 
five times that of the smaller platelet so the asymmetry of the mixture is large. The I-N 
spinodal remains close to the isotropic branch of the binodal; more so than in Fig. 7b, the 
(ci, c 2 ) representation of Onsager theory. The (x,p*) representation for FMT (Fig. 5c) shows 
that the general trend of the re-entrant bend moving to lower composition has continued, 
here reaching as low as about x = 0.2 for FMT. There is a very narrow region of coexistence 
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between the isotropic and nematic state of species 2 up to about p* = 0.06, suggesting that up 
to this pressure there is only a small fractionation effect. This fractionation becomes wider 
on increasing the pressure beyond this point, with the phase boundaries almost reaching 
x — 0, 1 by p* — 0.1. The re-entrant bend is also a dominating feature of the phase diagram 
in Fig. 7d for Onsager theory suggesting the FMT phase behaviour is correct, however it 
is not as pronounced as for FMT and the narrow I-N2 handle present at low pressures for 
FMT is wider for Onsager theory. 

In Fig. 8 we show the partial nematic order parameters S{ in the coexisting nematic 
phase(s) for A = 2,2.5,4 and 5, as obtained from FMT. For A = 2 (Fig. 8a) S\ and S2 
both rise rapidly near zero composition to near unity. S\ takes on smaller values than S 2 
for a given mole fraction as has already been observed for the lower size ratios. Hence the 
(smaller) particles of species 1 are much less ordered than particles of species 2 at a given 
composition, since the particles of species 2 are high in number and there is more freedom 
for the smaller particles to rotate in the dense system. Both partial order parameters remain 
close to unity as the composition x is increased from to about 0.5. The system is at N-N 
coexistence over this range (see Fig. 4c for the corresponding phase diagram). Stot reaches 
as low as about 0.5 and Si becomes as small as 0.2 at x — 1. For A = 2.5 (Fig. 8b) the 
N-N coexistence does not end at a critical point, leaving an interval (in composition) that 
the curves Si(x) do not enter. For A = 4 (Fig. 8c) S tot varies between about 0.3 and almost 
1. This range of Stot increases further for A = 5 (Fig. 8d) where it reaches as low as about 
0.2. For A = 4, 5, Stot & l so increases more rapidly as the composition increases than is the 
case for the lower size ratios. At high compositions, the (small) particles of species 1 possess 
very low nematic order, reaching less than 0.1 due to the larger size difference between the 
two species. 

Our phase diagrams share several features with those reported in previous studies of 
binary mixtures of anisometric hard core particles. Widening of the I-N phase coexistence 
region upon increasing the bidispersity parameter was found in binary mixtures of thick 
and thin rods [stJ, long and short rods 28], mixtures of rods and platelets 64], as well as 
binary mixtures of platelets usin g bo th the Parons-Lee scaling of the Onsager functional 
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651 ] and the Zwanzig model 36|, [37|. N-N phase coexistence ending in an upper critical 



point, as we find for an intermediate range of values of size ratios A, occurs in certain regimes 
for binary mixtures of rods and platelets. A notable difference to other studies of platelet 
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mixtures [65( is that we do not find N-N coexistence ending in a lower critical point, at 
least not for the range of densities and size ratios that we explored. Similar results to ours, 
where N-N coexistence ends in an upper critical point were obtained using Onsager theory 
for mixtures of thick and thin rods [27]. This system of rods also displays isotropic-isotropic 
phase coexistence at high enough values of the diameter ratio, which we do not find to be 
stable in the present study. 

IV. CONCLUSIONS AND OUTLOOK 

We have studied the bulk phase behaviour of binary mixtures of hard platelets, including 
isotropic and nematic states. The platelets are assumed to have circular shape and vanish- 
ing thickness. We have not considered positionally ordered phases such as columnar and 
crystalline phases, which are not expected to occur for the present model of particles with 
zero volume and zero packing fraction. Note that a first order N-Col transition was found 



at non-zero packing fraction in the limit of vanishing thickness 66[; simulations in this limit 



are possible because the model can be mapped onto a system of particles with finite volume 



but variable shape. Also, platelets with non-zero thickness exhibit a N-Col transition 



this was also reported for the Zwanzig model 35]. 



67j; 



For a variety of size ratios, we have compared the results of Onsager theory and FMT. 
In the monodisperse limit we find that Onsager theory overestimates the I-N transition 
densities and predicts a larger biphasic gap than FMT; the results for FMT in this limit 
compare quantitatively well with those from simulation. Both theories predict a first order 
I-N phase transition. We expect that Onsager theory overestimates the density jump at 
coexistence and also overestimates the size of the biphasic gap. The FMT results show the 
appearance of re-entrant phenomena at a lower size ratio than for Onsager theory. FMT 
results also show a larger range of compositions for which N-N demixing at a given size ratio 
occurs. The N-N demixing occurs at A = 2, where it is closed by an upper critical point. 
For A > 2.5 there is no critical point up to the densities we consider. In an experimental 
system of platelets with nonzero volume, one would expect a positionally ordered phase 
to be favoured before any possible remixing into a homogeneous nematic phase. We also 
examine the degree of nematic ordering along the phase boundaries for a selection of size 
ratios up to A = 2 where N-N demixing ends in an upper critical point. The partial nematic 
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order parameters S±, S2 vary smoothly with increasing composition and hence so does S to t- 
Where a re-entrant feature occurs, all three order parameters along the nematic branch of 
the binodal drop sharply. As A increases, S2 becomes smaller at x — 1, reaching < 0.2 for 
A = 2. We do not find any stable /-/ demixing for the range of A in this investigation. 

Through the range of A values considered, a striking observation is that whilst FMT 
predicts the occurence of phase boundaries at locations quantitatively different from Onsager 
theory, the topology of the phase diagrams for a particular choice of A is the same as for 
Onsager theory. Coupled with the fact that the I-N transition for the monodisperse case 
agrees well with simulation results |25[] , we gain confidence that the phase diagrams predicted 
by FMT are quantitatively more reliable since we assume that Onsager theory predicts the 
correct physics qualitatively. Furthermore, we have confidence that the order parameter 
profiles predicted by FMT are closer to those that would be observed from simulation studies, 
since the monodisperse limit in the present theory yields results for the I-N transition and 
nematic order parameter that is in good agreement with simulation results. Whether the 
differences to the true transition densities for the case of binary mixtures increases remains 
to be seen. 

The results from Onsager theory cannot be obtained by some scaling of the results that 
exist for binary mixtures of rods (thick and thin [27J or long and short [28]). For example, 
for there to be a mapping between thick and thin rods [23] and the present system, we would 
require simultaneously that (l+d) = A 2 +A and d = A 3 , where d = D2/D1 > 1 is the diameter 
ratio of the thick and thin rods. These requirements are obtained from consideration of the 
free energy for thick and thin rods (compare Eq. (2) of Ref. 23] with Eq. (fl3l) of the present 
investigation). Clearly these conditions cannot be satisfied simultaneously except for the 
trivial case A = d — 1. In the monodisperse limit, the mapping of concentrations from rods 
to platelets is c ro( j = (7r 2 /2)c p i ate iet where c ro d = (ti/4:)L 2 D is the conventional dimensionless 
concentration for rods (L is the rod length, D is the rod diameter). Clearly, simulation 
results for the systems studied in the present paper are most desirable. The predictions 
of the present work could also be tested experimentally. For example, for gibbsite (having 
typical experimental radius of approximately lOOnm) and hydrotalcite (typically between 
about 25 and 75 nm) Q] one has values of A between 1.3 and 4, which are in the range of 
the present work. We have taken mean radii here. Of course in experiments there will be 
effects due to polydispersity and our model ignores effects of finite platelet thickness. 
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Interesting future work could involve examining the phase behaviour of binary mixtures 
of polarizable platelets, that is particles that interact with some external applied magnetic 
field, as has been investigated for binary rod mixtures 
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of rods and platelets could be investigated 
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691 ] . Also, binary mixtures 



474|. Another line of investigation would 



be to incorporate polydispersity into the model since many platelet systems often have 

ead to an extension of the 



a significant polydispersity, see for example This would 
work by Speranza and Sollich and others on rod-like particles, 
Vroege on thickness-polydisperse platelets [781 ] . 



75M77j and by Wensink and 



Appendix A: Calculation of Relevant Integrals 

Here we evaluate some integrals that are important in the bifurcation analysis. The first, 
Ii, is a standard integral in excluded volume calculations 



Zi = J duj J do/ sin 7 = J duj J div'y/ 1 — (a? ■ a/) 2 . 



(Al) 



Without loss of generality, let uj' = e z in the integrand, where e z is the unit vector in the 
z-direction (the symmetry axis for the system). Now we have 



Xi = J duj' J duy/l- (w-e z ) 2 . (A2) 
Remembering that in our coordinate system uj' = (cos 9' sin 9', sin <p' sin 9', cos 9'), 



Xi = I duj' J cWl - cos 2 9 

p2tt P7V 

dJ \ d<p d9 sin 2 9 
Jo Jo 

2 / J, J _ /i„3 



= 7T 2 J dJ = 4tt 3 (A3) 

as required. The second integral arises in FMT. We wish to integrate the triple scalar 
product over three unit orientation vectors, uj,uj',uj". 



X 2 = J du J duj' J duj"\uj ■ {uj' x «")|. (A4) 
Similar to the X\ calculation, let uj" = e z in the integrand. 

J 2 = J duj" J duj J duj'\uj ■ («' x e z )\. (A5) 
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Remembering that in our coordinate system uj 
(cos 9' sin 9', sin (j)' sin 9', cos 9') , 



X 9 = / du" / / duj' 



( sinfl^ ' 





COS6 1 



sin 0, 0, cos 9) and a/ 



V 



y cos a j 
duj" I duj I du>'\ sin sin sin </>' 



sin <p' sin 9 
— cos 9' sin 0' 




dw" ■ 2tt 3 = 8tt 4 



(A6) 



as required. (Here the inner two integrals alone give 7r 3 /2 which can be found straight- 
forwardly). There are three further integrals required for the bifurcation analysis. X3 = 
J du> [P 2 (cos#)] 2 = 47r/5 trivially. Let us turn our attention to 



X 4 = J duj J dJ 'sin7[P 2 (cos0)P 2 (cos0')] 



We expand sin 7 as 



with coefficients 



C2n 



sin 7 = ^c 2 „P 2 „ (cos 7) 

n=0 

-yr(4n + l)(2n - 3)!!(2n - 1)!! 

2 n+1 (n + 1)! ' 



where the double factorial is defined by 

n- (n-2)---5-3- 1 n > odd 
"!! = < n- (n-2)---6-4-2 n > even 

1 71= -1,0. 

Therefore the inner integral of flA7j) becomes 



cfa/sin7P 2 (cos0') = / du' ^ c 2n P 2n (cos 7) P 2 (cos 9') 

J n=0 

/*7T 00 

d(f>' / tiff' sin 0' yj c 2n P 2n (cos 7) P 2 (cos fl'). 



We now utilise the addition formula for Legendre polynomials, 



2n 



P 2 „(cos7) = P 2n (cos0)P 2n (cos0') + 2 ^ 



(2n — rn)! 
(2n + to)! 

x P™(cosfl)P 2 ™(cosfl') cosm(0 - (j)'). 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



(A12) 
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Hence 



/PZ7T 
dd sin ^P 2 (cos 9') = 2tt J sin & d9' 

oo 

X ^ C 2n P 2n( COS 9) P 2n(<X>8 9')P 2 (COS 0) 



(A13) 



n=0 



where the sum involving associated Legendre functions of the first kind vanish in the inte- 
gration over 0'. Now we introduce the notation 



So 



(P 2n ) f = / d9'f(9')P 2n (co S 9'). 
Jo 



/oo 
cLi/sin7P 2 (cos6>') = 2n ^ c 2n P 2n (cos 0)(P 2n ) 
n=0 



f 



(A14) 



(A15) 



with / = P 2 (cos0'). From the orhogonality conditions of Legendre polynomials 79] we have 
that (P 2n )f = 2/(2m+ l)5 mn where 5 mn is the Kronecker-5 symbol. Here, m = 2 and, given 
that c 2 = — 57r/32, we have 

f 2 7T 2 

/ du/sin7P 2 (cos0') = 27rc 2 P 2 (cos0) • - = -— P 2 (cos0). 
J 5 2 

Hence 



(A16) 



j du)P 2 {co$6) 




-^P 2 (cos0) 







2tt 3 / d0sin0[P 2 (cosi 



-2tT 



1 7T 3 

5 ~10 



(A17) 



as required. The last integral to consider is 

X 5 = J d u J dd J du>"\u> ■ (d x d')\[P 2 (cos9)P 2 (cos9') 
+ P 2 (cos0)P 2 (cos0") + P 2 (cos0')P 2 (cos0")]. 

We consider integrating the first term of (1A18I) . namely 

X, = / duj [ dd [ dd'\w ■ (d x d')\P 2 ( y cos6)P 2 ( y cos9'). 



(A18) 



(A19) 



We let d' = e z in the integrand so that 

X 5 = / dd' j du: [ dd\ sin7(n ■ e z )|P 2 (cos 0)P 2 (cos 9') 



(A20) 
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where 7 is the angle between u> and u>' and n is a unit vector perpendicular to both uj and 
uj' . Since e z is fixed, n • e z — sin#". Hence 



X 5 comprises two other similar terms so by symmetry, X 5 = 3X 5 . Hence X 5 = — 37r 4 /5 as 
required. 
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FIG. 1: Results for A = 1.1. (a) shows the phase behaviour in the (x,c) representation, (b) in 

the (c\,C2) representation and (c) in the (x,p*) representation. The results for Onsager theory 

and FMT are plotted on the same graph. In (a)-(c) the upper pairs of solid lines are the binodals 

according to Onsager theory, the lower pairs of solid lines are the binodals according to FMT 

and the dashed lines indicate the I-N spinodals. An isotropic phase (J), isotropic and nematic 

coexistence (/ + N) and a nematic phase (TV) are present. In (c) the inset is for Onsager theory. 

In (a) and (b) the thiner solid lines are selected tie lines connecting coexisting state points. Note 

that the two theories are plotted on the same graph for comparative purposes; if one is interested 

in the FMT results for example, then the phase behaviour above the nematic binodal is nematic 

only and one should ignore the results from Onsager theory. In (d) we plot the nematic order 

parameters along the nematic binodal for Onsager theory and FMT. The dotted curve is Si, the 
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dashed curve is S2 and the solid line is <Stot- 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 3: Results for A = 1.7. The notation is the same as in Fig. 1. In (c) for FMT, if one keeps a 
constant fluid composition, for example at around x = 0.59 but increases the pressure from p* = 
to p* = 1 the state changes from I — > I+N2 — > N2 — > I+N2 — > N2. In (d) we present the order 
parameters only for FMT. 
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0.2 0.4 06 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1.2 




FIG. 4: Phase diagrams for A = 2. (a) shows the phase behaviour in the (ci,C2) representation for 
Onsager theory and (b) for FMT. (c) shows the phase behaviour in the (x,p*) representation for 
Onsager theory and (b) for FMT. The topology is the same for Onsager theory and FMT. Solid lines 
denote the binodals and the thinner solid lines denote tie lines. As well as I-N coexistence, there 
is now coexistence between two nematic phases; the N-N phase ends in a critical point (depicted 
as a shaded circle) in both Onsager theory and FMT, though the critical point is at lower densities 
for FMT. The dashed lines denote the I-N spinodals but we do not show the location of the N-N 
spinodals on these phase diagrams. There is also I-N-N coexistence. In (a) and (b) this triple 
point is represented by a region bounded by dotted lines, the vertices of which are the pure phases. 
In (c) and (d) this region is collapsed onto a triple line, shown here with a dotted line. 
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x x 



FIG. 5: Phase diagrams for A = 2.5. The notation is the same as in Fig. SJ The N-N coexistence 
does not end in a critical point in these phase diagrams. 
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FIG. 6: Phase diagrams for A = 4. The notation is the same as in Fig. 01 We have focused on 
the lower part of the phase diagram as this is where the most interesting phase behaviour occurs. 
The N-N phase separation continues to high densities beyond the scale of these plots, with a huge 
immiscibility gap. 
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FIG. 7: Phase diagrams for A = 5. The notation is the same as in Fig. [31 We focus on the re-entrant 
behaviour in this figure. The general trend of the re-entrant bend moving to lower composition 
continues here once again. 
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FIG. 8: The nematic order parameters along the nematic binodal for A = 2,2.5,4 and 5 using 
FMT. The dotted curves indicate Si, the dashed curves indicate 52 and the solid curves indicate 

Stot- 
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